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1.  INTRODUCTION 

In  a  tropical  cyclone,  the  concentrated  potential  vor- 
ticity  (PV)  source  due  to  heating  in  a  circular  eyewall  can 
result  in  a  reversal  of  the  radial  PV  gradient,  allowing  the 
vortex  to  become  barotropically  unstable.  In  this  manner  an 
axisymmetric  vortex  can  develop  asymmetries,  redistribute 
PV  through  chaotic  nonlinear  mixing,  and  eventually  resym¬ 
metrize  with  a  different,  stable  structure.  Recent  studies 
of  this  process  within  the  framework  of  an  unforced  non- 
divergent  barotropic  model  (Schubert  et  al.,  1998,  herein 
referred  to  as  S98)  have  helped  provide  insight  into  diverse 
aspects  of  hurricane  dynamics,  including  the  development 
of  spiral  bands  and  mesoscale  vortices,  the  existence  of 
polygonal  eyewalls,  and  asymmetric  eye  contraction. 

One  approach  to  understanding  the  PV  redistribution 
process  is  to  treat  it  as  an  initial  value  problem.  Numeri¬ 
cal  results  from  nonlinear  models  can  simulate  the  process 
for  a  limited  time  but  cannot  track  the  details  indefinitely, 
since  the  PV  field  tends  to  filament,  developing  structure 
on  increasingly  finer  scales  that  eventually  cannot  be  re¬ 
solved.  An  alternative  approach  is  to  predict  directly  the 
final  steady  flow  produced  by  the  mixing,  without  simu¬ 
lating  the  time  evolution  of  the  flow.  At  least  two  ap¬ 
proaches  have  been  proposed,  based  on  minimum  enstro- 
phy  and  maximum  entropy  arguments.  This  paper  concen¬ 
trates  on  the  minimum  enstrophy  vortex  (MinEV)  problem, 
generalizing  previous  results  for  the  unforced  nondivergent 
barotropic  model  and  extending  the  analysis  to  the  diver¬ 
gent  barotropic  (shallow  water)  model. 

2.  NONDIVERGENT  BAROTROPIC  MODEL 

In  the  unforced  nondivergent  barotropic  model  the  en¬ 
strophy,  energy,  and  angular  momentum  are  all  conserved 
in  inviscid  axisymmetric  flow.  With  dissipation,  none  of  the 
three  are  conserved,  but  the  energy  and  angular  momen¬ 
tum  typically  decay  much  more  slowly  than  the  enstrophy. 
Thus,  at  least  for  intermediate  (non-diffusive)  time  scales, 
the  flow  may  be  thought  to  evolve  toward  a  state  of  min¬ 
imum  enstrophy  while  conserving  the  energy  and  angular 
momentum.  Starting  from  this  “selective  decay  hypothesis" 
it  is  possible  to  derive  the  final  state  of  the  flow  by  vari¬ 
ational  techniques,  yielding  the  MinEV  solutions  of  Leith 
(1984)  and  S98.  In  this  section  we  extend  the  approach 
of  S98  to  include  two  mixing  radii  and  to  constrain  both 
energy  and  angular  momentum  in  the  same  analysis. 
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Consider  an  axisymmetric  initial  flow  with  tangential 
velocity  vo(r)  and  corresponding  vorticity  Co  =  d(rvo)/rdr. 
If  the  flow  is  barotropically  unstable,  it  may  evolve  through 
chaotic  nonlinear  mixing  and  asymptotically  approach  a 
steady  final  state  which  is  again  axisymmetric.  We  assume 
that  the  mixing  is  confined  to  the  region  0  <  a  <  r  <  b, 
and  the  flow  outside  this  region  is  undisturbed.  Denoting 
the  tangential  velocity  and  vorticity  of  the  final  state  by 
u(r)  and  £  =  d(rv)/rdr ,  respectively,  the  enstrophy  excess 
of  the  final  state  is 

Z[v,a,b]  :=  f  5  (C2  -  Co)  rdr,  (1) 

J  a 

where  we  have  omitted  the  multiplicative  constant  27t  for 
convenience.  The  corresponding  energy  and  angular  mo¬ 
mentum  excesses  are 


E[v,a,b]  :=  f  | 

J  a 

( v 2  —  Vo )  rdr 

(2) 

and 

pb 

A[v,a,b]:=  /  r 

(v  —  vo)  rdr . 

(3) 

We  seek  v(r),  a,  and  b  which  minimize  Z[v}  a, 
the  constraints 

6]  subject  to 

E[v ,  a,  b]  =  0, 

A[v ,  a,  b]  =  0 

(4) 

and  the  boundary  conditions 

v(a)  =  v0(a)5 

v(b)  =  vo(b). 

(5) 

Using  standard  techniques  from  the  calculus  of  vari¬ 
ations,  we  find  that  the  first  variations  of  the  boundary 
conditions  (5)  are 


K(a)  ~  Co(a)]  Sa  +  Sv(a)  =  0.  (6) 

and 

[<(b)-to(b)]6b  +  6v(b)  =  Q.  (7) 

Using  these,  we  can  write  the  first  variation  of  Z  as 

6Z  =  -J'f5v rdr  ~  5  b5b  (8) 

+  |  [C(a)  —  Co(a)]2  a5a. 

The  first  variations  of  E  and  A  are  found  in  a  similar  manner 
to  be 

fb 

SE  =  /  vSvrdr  (9) 
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and 


3.  DIVERGENT  BAROTROPIC  MODEL 


6A  =  rSvrdr ,  (10) 

J  a 

where  in  both  cases  the  boundary  terms  vanish  due  to  (5). 

To  derive  equations  for  the  MinEV  solution,  we  note 
that  if  v ,  a,  and  b  minimize  Z  subject  to  the  constraints  of 
constant  E  and  A,  then  they  also  minimize  the  functional 

I[v ,  a,  b]  :=  Z[vy  a,  b]  -f  f3E[v ,  a,  b]  +  ^fA[vy  a,  6],  (11) 


known  as  the  Lagrangian,  where  j3  and  7  are  Lagrange 
multipliers.  We  can  combine  (8)-(10)  with  (11)  to  obtain 

61  =  (— —■  +  Qv  +  77^  Svrdr  ,  x 

Ja  \  dr  J  (12) 

-  |  K(6)  -  Co (b)f  bdb  +  I  K(a)  -  Co(a)]2  a<5a. 

Since  any  stationary  solution  must  have  6 1  =  0  for  all  vari¬ 
ations,  we  can  obtain  equations  for  the  stationary  solution 
as  follows.  First,  consider  variations  for  which  6a  =  0  and 
6b  =  0  but  6v  is  otherwise  arbitrary.  Then  61  =  0  im¬ 
plies  that  the  integrand  in  (12)  must  vanish,  leading  to  the 
Euler-Lagrange  equation 


d  d(rv) 
dr  rdr 


—  j3v  =  7  r, 


a  <  r  <b, 


(13) 


where  we  have  substituted  for  £  in  terms  of  v.  The  solu¬ 
tion  of  (13)  subject  to  (5)  may  be  obtained  analytically  in 
terms  of  Bessel  functions  and  the  parameters  a,  b,  and 
7.  To  determine  these  parameters,  we  substitute  (13)  into 
(12)  and  consider  cases  where  6a  =  0  and  6b  =  0  indepen¬ 
dently,  concluding  that  the  condition  61  =  0  also  implies 
the  transversality  conditions 


C(o)  =  Co  (a),  C(b)  =  Co  (b).  (14) 


The  conditions  (14)  and  the  constraints  (4),  i.e., 


f  \  {v2  ~  vo)  rdr  =  0, 
da 


(15) 


and 


v  —  vo)  rdr  =  0, 


(16) 


then  give  four  conditions  on  a,  6,  f3t  and  7.  The  cases 
treated  in  S98  are  recovered  in  the  limit  as  a  ->  0  with 
(5  =  0  (MinEV-M)  and  7  =  0  (MinEV-E). 

To  see  whether  the  stationary  solution  defined  by  (13)- 
(16)  in  fact  corresponds  to  minimum  enstrophy  (rather  than 
a  maximum  or  saddle  point),  we  compute  the  second  vari¬ 
ation  of  I  (at  the  stationary  solution),  obtaining 


<s2/  =  i 


d(6v) 


dr 


(l3+±)(6v)2 


rdr .  (17) 


From  (17)  we  conclude  that  for  nontrivial  continuous  vari¬ 
ations  6v  satisfying  6v(a)  =  6v(b)  =  0  we  have  62I  >  0, 
provided  that 

P>~^-  (18) 

Thus — if  the  condition  (18)  holds — the  stationary  solution 
defined  by  (13)— (14)  and  (4)  corresponds  to  a  minimum  of 
I  and  thus  is  indeed  the  minimum  enstrophy  vortex. 


The  MinEV  analysis  for  the  axisymmetric  divergent 
barotropic  (shallow  water)  model  on  an  /-plane  parallels 
that  developed  above  for  the  nondivergent  model.  Since 
the  free  surface  height  h  is  variable  in  this  model,  the  key 
quantity  conserved  in  inviscid  flow  is  the  potential  vorticity 


P  = 


H 

h 


f  |  d(rv) 
J  rdr 


(19) 


where  H  is  a  constant  reference  height.  We  assume  that 
the  wind  and  mass  fields  are  always  in  gradient  balance, 
i.e., 


where  g  is  the  acceleration  due  to  gravity.  Since  the  po¬ 
tential  enstrophy,  mass,  total  energy,  and  absolute  angu¬ 
lar  momentum  are  all  conserved  in  inviscid  flow,  the  se¬ 
lective  decay  hypothesis  suggests  that  the  flow  may  evolve 
toward  a  state  of  minimum  potential  enstrophy  while  ap¬ 
proximately  conserving  the  mass,  absolute  angular  momen¬ 
tum,  and  total  energy.  The  problem  is  in  one  sense  simpler 
than  the  nondivergent  barotropic  case,  since  the  global  de¬ 
pendence  of  the  mass  and  wind  field  on  the  potential  vor¬ 
ticity  (through  the  inverti bility  principle)  implies  that  there 
is  no  “mixing  radius" :  the  problem  should  be  treated  on 
the  whole  domain  0  <  r  <  00. 

We  can  write  the  potential  enstrophy,  mass,  total  en¬ 
ergy,  and  angular  momentum  in  deviation  form  as 


r  00 

Z[h,v]:=  |(  P-ffhrdr,  (21) 

^0 


poo 

M[h]  :=  I  (h  —  H)  rdry 
Jo 

POO 

E[h ,  v]  :=  1  [v2h  +  g(h  -  H)2]  rdr , 

Jo 


(22) 


(23) 


and 


POO 

A[h,v]:=  /  [rvh  -h  \fr2(h  —  77)]  rdr.  (24) 

Jo 


We  assume  that  P  =  /  outside  a  bounded  region  (and  that 
h  is  bounded);  then  it  can  be  shown  that  v  and  h  —  H  both 
decay  like  0(e“r)  as  r  oo,  so  each  of  Z,  M,  E,  and  A 
are  finite.  Thus,  the  MinEV  problem  is  to  find  the  flow  field 
h,  v  in  gradient  balance  which  minimizes  Z[hy  v]  subject  to 
the  constraints  M[h]  =  M[h0],  E[hyv]  =  E[h0,v 0],  and 
A[hy  v]  =  A[hoy  uo],  where  ho  and  vo  specify  the  initial  flow 
(assumed  to  be  axisymmetric  and  in  gradient  balance). 
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Using  standard  techniques  from  the  calculus  of  varia¬ 
tions  we  can  compute  the  first  variations  of  Z[h,v],  M[h], 
E\h,v],  and  A[h:  v]  using  (19)  and  (21)-(24),  obtaining 

sz  =  fo  [~Hp2  -  f2)Sh~  H<^:5v]  rdr’  (25) 


r  oo 

SM  =  /  Shrdr) 
Jo 


(26) 


POO 

SE=  {[|v2  +g(h-  H)]5h  +  vh6v}  rdr,  (27) 
Jo 

d 

r  oo 

SA  =  ({rv  A  \}r2)Sh  A  rhSv\  rdr.  (28) 
Jo 


To  minimize  Z  while  holding  M,  E ,  and  A  constant,  we 
introduce  the  Lagrange  multipliers  a,  0,  and  7  and  seek  to 
minimize  the  functional 


I[ht  v]  :=  Z[h,  v]  +  aM[h]  A  0E[h ,  v]  +  7 A[h,  v].  (29) 
Using  (25)-(28)  in  (29)  we  obtain  the  first  variation 

SI  =  J  [-UP2-f2)  +  a  +  p{\v2  +  g{h-H)) 

A  7  (rv  A  ^fr2)  J  Sh  rdr 

r  n 

I  —H  — - b  0vh  A  77-/1 1  Sv  rdr , 

(30) 

which  must  vanish  if  I  is  stationary  at  the  solution  h,  v. 

Before  proceeding  further,  we  note  that  Sv  and  Sh  are 
not  independent  variations,  since  they  are  related  by  the 
constraint  of  gradient  balance.  We  can  construct  a  corre¬ 
sponding  functional  by  multiplying  (20)  by  an  unspecified 
function  p(r)  and  integrating,  obtaining 

G[h,v]  :=  f  [{f+l)v~9^:]  Vrdr. 

J  0 


Since  we  must  have  SI  —  0  for  the  stationary  solution 
h,  v ,  we  can  obtain  Euler-Lagrange  equations  from  (32)  as 
follows.  First,  since  the  function  fi(r)  was  arbitrary,  we  can 
choose  it  so  that  the  term  in  brackets  in  the  first  integral 
in  (32)  vanishes,  i.e., 

(33) 

-7  (rv  A  \  fr 2)  . 

Then  the  first  integral  in  (32)  vanishes,  regardless  of  Sh. 
Thus,  the  second  integral  must  also  vanish  for  all  variations 
Sv,  so  we  conclude  that  the  term  in  brackets  in  that  integral 
must  be  zero,  i.e., 

=  Pvh  +  rrh+  (f  +  y)v>  (34) 

These  two  equations,  along  with  the  definition  of  P 

Hd{rv)  =hp_  Hf  (35) 

rdr 

and  the  gradient  wind  equation 


form  a  system  of  ordinary  differential  equations  to  be  solved 
for  the  four  variables  /x,  P,  v,  and  h.  The  boundary  condi¬ 
tions  on  this  system  are 

t»(0)  =  o,  MO)  =  0  (37) 

(which  come  from  examining  the  behavior  of  Taylor  approx¬ 
imations  to  the  solution  as  r  — >  0)  and 

lim  P(r)  =  /,  lim  h(r)  =  H.  (38) 

r-4oo  r— ►oo 

Finally,  the  constraints  on  M,  E ,  and  A  serve  to  determine 
the  values  of  the  Lagrange  multipliers  a ,  0 ,  and  7  in  terms 
of  the  initial  mass  and  velocity  ho  and  Vo. 
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■jOt) 
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rdr 
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